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ABSTRACT 



Although the 55 Cnc system contains multiple, closely packed planets that are 
presumably in a coplanar configuration, we use numerical simulations to demon- 
strate that they are likely to be highly inclined to their parent star's spin axis. 
Due to perturbations from its distant binary companion, this planetary system 
precesses like a rigid body about its parent star. Consequently, the parent star's 
spin axis and the planetary orbit normal likely diverged long ago. Because only 
the projected separation of the binary is known, we study this effect statisitically, 
assuming an isotropic distribution for wide binary orbits. We find that the most 
likely projected spin-orbit angle is ~50°, with a ~30% chance of a retrograde 
configuration. Transit observations of the innermost planet - 55 Cnc e - may 
be used to verify these findings via the Rossiter-McLaughlin effect. 55 Cancri 
may thus represent a new class of planetary systems with well-ordered, coplanar 
orbits that are inclined with respect to the stellar equator. 



1. Introduction 



The 55 Cancri planetary system contains five known planets in low-eccentricity or- 
bits between 0.01 and 5.75 AU ( Butler et aDll997l : iMarcv et al]l2002[ iMcArthur et al.lbood : 



Fischer et al.ll2008l ). Recently, the inne r most planet in the s ystem - the ~8 55 Cnc e 



was found to transit ( IWinn et al.ll201ll : iDemory et al.ll201ll ). This is particularly intriguing 
because it may allow observers to measure (among other pr operties) the spin-orb it angle of 
the planetary system via the Rossiter-McLauglin effect (see lCaudi &: WinnI 120071 ). 
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Rossiter-McLauglin measurements have revealed that the orbits of many hot Jupiters 
are hi ghly inchned relative to their host star's spin axis and are eve n retrograde in some cases 



ir s sp 
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[e.g., [Anderson et al.ll2UlUI : liriaud et al.ll2UlUt iWinn et al.ll2010l ). This conflicts with stan- 
dard models of planet migration (ILin et al.lll996[ ) because protoplanetary disks are generally 



though t to be aligned with the stellar equator, as has been confirmed by debris disk measure- 
ments (jWatson et al.ll201ll ). The origin of highly inclined or retrograde planetary orbits may 
therefore be dynamical. Several recent papers have shown that such orbits may be produced 
by tida l circularization of high ly eccentric/inclined planets after a phase of planet-planet scat- 
tering ([Nagasawa et al. 2008 ) and/or Kozai interactions with a binary star or massive pl anet 



( iFabrycky &: Tremaind 12007 



Naoz et al.ll201ll : iKatz et al.ll201ll : iLithwick fc Naozll2011f ). 



Another noteworthy feature of 5 5 Cnc is that it also h as a 0.27 Mq binary companion 



at a projected distance of 1065 AU (iMugrauer et ahlbooeh. While there have been many 



dyna r nical studies of this planetary systern (e.g., [Raymond et al.ll2006l : iBarnes fc Greenberg 
20071 : iFischer et al.l l2008l : [Raymond et al.ll2008l ). none have investigated the effects of the 
binary companion. Presumably, this is because this well-ordered system displays none of the 
excited orbits expected from severe impulses or Kozai interactions with its binary companion. 
Here we show that this binary companion does in fact play a major role in the planetary 
system's dynamics. Using numerical simulations, we demonstrate that it is very likely this 
binary has driven the planetary system to a very high (and possibly retrograde) spin-orbit 
angle with respect to the parent star's spin axis. Moreover, the Kozai oscillations studied in 
recent works do not produce this this high spin-orbit angle. Rather, it is a mechanism first 
described in llnnanen et al.l (119971 ). where the entire planetary system smoothly precesses as 
a rigid body. 



Numerical Methods 



To model the dynamics of the 55 Cancri system, we use the wide binar y algorithm 
of the MERCURY integration package (IChamberd Il999l : [Chambers et al.l I2OO2I ) . All of our 
simulations include both the primary and secondary stars of 55 Cnc as well as one or more 
of the syst em's known planets. In each simulation, the primary's mass (55 Cnc A) is set to 
0.90 5 Mr7^ (Ivon Braun et al.ll201ll ). while the secondary's mass (55 Cnc B) is fixed at 0.27 
Mq (IMugrauer et al.ll2006l ). In total, we perform 501 simulations. 450 of our runs include 
55 Cnc d, the most distant and massive planet in the system, as well as a hypothetical 
Saturn-mass planet we designate as 55 Cnc g. Our motivation for using such systems is 
discussed in Section 3. These simulatio ns are integrated for 1 0.2 Gyrs, as this matches the 
most recent age estimate of the system (Ivon Braun et al.ll201ll ). These runs are divided into 



- 3 - 



nine 50-simulation subsets. In each subset, 55 Cnc B is assigned a different semimajor axis: 
750, 1000, 1250, 1500, 2000, 3000, 4000, 5000 or 8000 AU. All of its other orbital elements 
are drawn randomly from an isotropic distribution. 

In our remaining simulations, we repeat some of our 2-planet simulations using a differ- 
ent planetary configuration. For 48 of these reruns, we use the four outermost planets (55 
Cnc b, c, f, and d) in place of 55 Cnc g and d. These simulations are only run for 1 Gyr 
because the orbit of 55 Cnc b necessitates a timestep of just 1 day. Finally, we also rerun 
3 of these simulations using all 5 planets. These 5-planet simulations are run for only 50 
Myrs, as they require a timestep of ~1 hour. The initial orbits and masses assumed for 55 
Cue's planets are shown in Table 1. 

In addition to the planets' and stars' mutual gravitational interactions, we also include 
external perturbations from the Galactic tide and passing field stars, since these forces affect 
the dynamics of 55 Cnc B for s emimajor axes beyon d a few thousand AU. We employ the 
Galactic tidal field described in lLevison et al.l (120011) . To m odel perturbations from passing 
stars, we use the impulse approximation (IRickmanl Il976l) . These impluses are generated 
by randomly drawing from a stellar mass functio n (IReid et al.l 120021 ) and a mass-dependent 
velocity distribution ( iGarcia-Sanchez et al.ll200ll ). 



3. Results and Discussion 



In Figure 1, we display the inclination evolution seen in one of our simulations. This 
simulation is one of our 4-planet simulations, and 55 Cnc B is placed on an initial orbit 
with a = 1250 AU, e = 0.93, and i = 115°. (Based p rojected separation, the most probable 
semimajor axis is ~1340 AU (IFischer fc Marcylll992l ).) We see in this figure that the planets 
of 55 Cnc leave their initial orbital plane very rapidly. In fact, within 30 Myrs the planets 
have reached a retrograde configuration relative to their original orbital plane. If the parent 
star's spin axis is initially aligned with the planets' orbits, this means that the spin-orbit 
angle of the system quickly diverges from 0°. For this entire Gyr simulation, the planets 
continue to precess between prograde and retrograde configurations in a very regular manner. 
Interestingly, the inset plot of this figure illustrates that all four planets maintain a very tight 
coplanar configuration throughout this evolution. 



m 



The dynami c al me chanism that drives this inclination evolution was first documented 



Innanen et al.l (Il997( ). In a system consisting of just one planet embedded within a binary 



star system, perturbations from the secondary star will cause the planetary orbit to undergo 
Kozai-like behavior, where the inclination and eccentricity of the planet oscillate exactly out 
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of phase as the longitude of pericenter (cu) circulates or librates ( lKozailll962[ ). The frequency 
of this oscillation is a function of the star's mass, semimajor axis, and eccentricity as well 
as the planet's semimajor axis. However, when more planets are added to the system the 
behavior can change. This is because the evolution of u may no longer be dominated by the 
stellar companion's perturbations. Instead, the mutual interactions between the planets can 
drive the evolution of the longitude of pericenter. As long as the precession timescale of u 
is much shorter than the Kozai timescale due to the binary's perturbations, the evolution 
of the planetary orbits will not resemble a Kozai resonance. Indeed, integrations of the 55 
Cnc planets in isolation indicate that u of planet d circulates every 2.5 Myrs. In contrast, 
integrations with only planet d and a binary companion (a = 1340 AU, e = 0.95) yield a 
much larger typical circulation period of 10^~® Myrs. Thus, the Kozai mechanism will not 
operate, and the planets' eccentricities and mutual inclinations remain low. 

However, in addition to a;, the lo ngitude of ascending node also precesses in the 
reference frame of binary orbital plane (jinnanen et al.lll997l ). Unless the binary star's orbital 
plane and the initial planetary orbital plane coincide, the precession of in the binary's 
reference frame will translate to an inclination precession with respect to the initial planetary 
plane. Although each planet would have a different Vt precession rate i n isolation, the self- 



gravity of the syste m causes the planets to precess at a uniform rate (ITakeda et al.l 12008 



Batygin et al.l l201ll ). Consequently, the system maintains a rigid, coplanar shape, even 
though its orbital plane can become greatly inclined to the star's spin axis (assuming this 
axis is perpendicular to the original planetary orbital plane). 

Figure 1 also shows that although our simulated system spends time in retrograde con- 
figurations, its inclination never exceeds a certain value, imax (in this case imax — 130°). The 
reason for this is that the vertical component of the planetary orbital angular momentum, 

is always observed to be conserved in the reference frame of the binary orbital plane. 
Another way of saying this is that as precesses, the angle between the orbit normal of 
the planetary system and the orbit normal of the binary is fixed. Thus, the inclination with 
respect to the original planetary orbital plane reaches a maximum after VL has precessed 180° 
in the binary orbital frame. In the reference frame of the original planetary orbital plane 



COSimax = COs2iun (1) 

where ibin is the inclination of the binary with respect to the planetary orbital plane. 

In the above simulation, we do not include the innermost planet of the 55 Cnc system, 
planet e. Because of its tiny orbital period, the computing costs to follow such an integration 
for 1 Gyr would be prohibitive. However, this planet is only 0.1 AU from the next nearest 
planet (b), and it is subject to powerful perturbations from interplanetary gravitational 
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interactions just as all the other planets are. Consequently, it too should conform to the 
rigid body nature of the evolution displayed in Figure 1. 

Fortunately, verifying this does not actually require a 1-Gyr integration. A few of our 
4-planet simulations have inclination precession rates even faster than the system shown 
in Figure 1 (due to a lower pericenter and semimajor axis for the binary orbit). These 
configurations with faster precession rates allow us to use shorter numerical simulations to 
verify that planet e follows this rigid body precession. In Figure 2, we display one such 
case. This figure actually displays inclinations from two different simulations: planet e's 
inclination in a 5-planet simulation and the inclinations of planets b, c, f, and d from a 
4-planet simulation. The same binary orbit is used in both simulations. We see that the 
behavior of each system is very similar, with inclinations oscillating between 0° and ~95°. In 
addition, the period of oscillation is nearly the same, although the 5-planet system precesses 
a little slower (due to the fact that planet e has a smaller "natural" Q precession rate and 
slows the mean rate down slightly). Because planet e is so close to its parent star, we include 
general relativistic precession as well as a J2 component to the parent star's potential ( J2 = 
5 X 10"''). However, the behavior of planet e does not change noticeably whether we include 
these effects or not. Thus, we conclude that the inclination behavior seen in our 4-planet 
simulations is an excellent proxy for planet e's inclination as well. 

Although we can use our 4-planet simulations to extrapolate the behavior of all 5 planets 
in 55 Cnc, this is still not an ideal computing situation. Because of planet b's short orbital 
period, even our 4-planet simulations use a timestep of ~1 day, and we would like to integrate 
our planetary systems for 10 Gyrs rather than 1 Gyr. In addition, we would like to evolve 
many hundreds of systems, each with a different binary star orbit, since only the binary's 
projected separation is known. 

We now argue that 55 Cue's inclination evolution can be modeled accurately using 
integrations that include just the outermost planet (d) accompanied by a fictitious inner 
planet (g), rather than the known configuration. The reason for this is that planet d contains 
most of the total planetary mass of 55 Cnc. Therefore, planet d dominates the system's self- 
gravity that maintains its rigid behavior. The other planets should almost behave as test 
particles being dragged along with planet d. The only role that the other planets play in 
this evolution is perturbing planet d's longitude of pericenter to prevent Kozai oscillations. 
Consequently, the inclinations of any configuration of inner planets will evolve similarly as 
long as they are stable and perturb planet d's oj sufficiently. Considering this, we choose 
to replace the inner 4 planets with a S aturn- mass planet at 3.5 AU (which we call planet 



g). Such a planet is known to be stable ([Raymond et al.ll2008l ). and its large semimajor axis 



enables us to increase our integration step from 1 day to 125 days. 



- 6 - 



To demonstrate the accuracy of these new simulations, we compare our 48 4-planet 
simulations with simulations using planet g in place of b, c, and f. For both our 4-planet 
simulations and our 2-planet simulations, we measure the maximum planetary inclination 
attained during the first Gyr. These are plotted against each other in Figure 3a. We see 
that there is a tight 1:1 correlation between the inclinations attained in both simulation 
sets. Another parameter we can compare between the 2-planet and 4-planet systems is the 
precession rates of their planetary inclinations. This is done in Figure 3b where we plot the 
precession periods (obtained from an FFT) from both simulation sets. We see that in general 
the two precession periods are very near each other, although the tighter configuration of 
the 4-planet system does systemically yield a slightly lower precession rate. Given that the 
inclination behaviors of these two simulation sets are nearly identical, we conclude that our 
2-planet simulations are a suitable model for the inclination of the real 55 Cnc system. 

In total we perform 450 2-planet simulations that are integrated for 10.2 Gyrs. Each 
simulation includes a binary companion set on a different randomly generated orbit. In a 
small fraction of our simulations (16%), binary perturbations destabihze the planets, ejecting 
one or both. Because this has not occurred in the real 55 Cnc system, we ignore these runs. 
Using only the stable runs, we sample the inclination of planet d relative to its initial plane 
every Myr between 7 and 10.2 Gyrs. The median values of those inclination samplings are 
plotted as a function of binary semimajor axis in the upper panel of Figure 4. In addition, 
we use error bars to mark the boundaries of the upper 10% and bottom 10% of inclination 
measurements in each semimajor axis bin. 

If we assume that the parent star's spin was originally aligned with the planetary orbits, 
this plot shows the current spin-orbit angle of 55 Cnc. We see that for binary semimajor 
axes below 4000 AU the most likely spin-orbit angle should be ~65°. Beyond semimajor 
axes of ~5000 AU, the precession timescale becomes much larger than the system's age, 
and the influence of the binary star wanes. Hence, for the very largest binary separations, 
the planetary system is likely to be found only at low inclinations. Lastly, we note that 
although our median inclination is always prograde, 29% of the inclinations recorded in our 
simulations are retrograde for binaries with a < 5000 AU. 

In the bottom panel of Figure 4, we show the cumulative probability distribution for 
the possible semimajor axis of 55 Cnc B. To calculate this distributi on, we assumed tha t 



the distribution of wide binary semimajor axes is uniform in log space (jPoveda et al.l 120071 ). 
Furthermore, we assumed that for a given semimajor axis, the distribution of all other orbital 
elements for wide binaries is isotropic. Based on this orbital distribution, we measured the 
relative fraction of time that projected separations of 1065 AU occur. We have seen in the 
upper panel that the binary companion has its strongest effects on the spin-orbit angle for 
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a < 5000 AU. Our probability distribution in the bottom panel indicates that there is a 
~95% probability that the semimajor axis of 55 Cnc B is within this range. Thus, it is very 
likely that 55 Cnc B has significantly altered the spin-orbit angle of this system. 

These simulations indicate that the current spin-orbit angle of 55 Cnc should be quite 
high and perhaps even retrograde. This result is particularly exciting because 55 Cnc e was 
just recently discovered to be transiting (jWinn et al.ll201ll : iDemory et al.ll201ll ). During this 
planet's 95-minute transit, we estimate a 48 cm/s magnitude for the Rossiter-McLaughlin 
effect, assuming the s tar's maximum rotational velocity is damped by the planet-to-stellar 



radius ratio squared (Ivon Braun et a. 



current detection limits (IMayor et al 



20111 : iGaudi fc Winn] l2007f ). This is slightly above 
201ll ). and it should be possible to measure the spin- 



orbit angle. Because the measured angle is only projected, we also display the range of the 
projected angle (A) in the upper panel Figure 4, assuming a random orientation on the sky. 
Although the median values of A are lower than the true i nclina tions, they are still well above 
the 30° "misaligned threshold" specified in lTriaud et al.l (120101 ) for all binary semimajor axes 
besides a = 8000 AU. 



4. Conclusions 



We show that it is very likely that perturbations from 55 Cnc B have significantly altered 
the spin-orbit angle of the planetary system of 55 Cnc A, even though the planetary system's 
self-gravity has preserved the planets' initial coplanarity. Hence, even very distant binary 
companions such as 55 Cnc B can substantially alter planetary architectures. Assuming wide 
binary orbits refiect an isotropic distribution, we demonstrate in Figure 4 that the most likely 
value of 55 Cue's projected spin-orbit angle is ~45-55°. Furthermore, there is a significant 
(29%) chance of a retrograde configuration for 55 Cnc. Because 55 Cnc is a closely packed, 
coplanar system, it would be unique among known highly inclined planetary systems. This 
feature could distinguish retrograde planets produced from rigid body precession vs. those 
generated from Kozai interactions and/or scattering events. Since planet e of this system 
transits its parent star, it may be possible to soon verify our findings using the Rossiter- 
McLauglin effect. 
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Table 1: List of planet masses, semimajor axes, eccentrici ties, and inclinations used in our 
simulations' initial conditions ( iDawson fc FabryckyI 120101 ). All other unlisted orbital ele- 
ments were randomly drawn from a uniform distribution from to 27r. Note that planet g 
is a stable fictitious planet added to some of our runs. 



^Although iDawson fc Fabrvckv (|2010l) list the eccentricity of planet e as 0.06, their analysis is also 
consistent with zero eccentricity, and it is likely such a closely orbiting planet would be tidally circularized 
rapidly. We have run a subset of alternative simulations with an eccentric planet e, and our results are 
independent of this orbital parameter. 

^To avoid 2-D integrations, we give each planet a slightly different random initial inclination between 
and 1 degree. 
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Fig. 1. — Evolution of the inclination of the outer four planets of the 55 Cancri system vs. 
time. Inclination is measured relative to the initial planetary orbital plane. The inset plot 
resolves the evolution each planet's inclination. 
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Fig. 2. — Inclination evolution of the 55 Cancri planets in two different simulations using the 
same binary orbit. The solid line shows the evolution of planet e in a simulation including 
all 5 planets. The data points plot the inclinations of the outer four planets in a simulation 
including just those planets. Inclinations are measured relative to the initial planetary orbital 
plane. 
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Fig. 3. — a: Maximum planetary inclination measured in 4-planet simulations vs. the 
maximum planetary inclination measured in 2-planet simulations. Inclinations are measured 
relative to the initial planetary orbital plane, b: The inclination precession period measured 
in 2-planet simulations vs. the precession period measured in 4-planet simulations. Data for 
both plots only come from the first Gyr of each simulation. 
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Fig. 4. — Upper panel: The predicted range of the spin-orbit angle of the 55 Cnc planetary 
system as a function of the system's binary semimajor axis. The circular data points mark 
the median value of the true spin-orbit angle, and the triangular data points mark the median 
value of the projected angle (A) that could be measured via the R-M effect. The error bars 
mark the range between the lower 10% and upper 90% of all possible angle values seen in our 
simulations. The ranges of the projected angle and true angle are offset slightly at each binary 
semimajor axis to avoid confusion. Bottom panel: Cumulative probability distribution of 
the possible semimajor axis of 55 Cnc B if its projected separation distance is 1065 AU. This 
distribution assumes that wide binary semimajor axes are uniformly distributed in logo- 
space. All other binary orbital elements are assumed to reflect an isotropic distribution. 
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